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Abstract 

An approach to the model-independent searching for the Z' gauge 
boson as a virtual state in scattering processes is developed. It ac- 
counts for as a basic requirement the renormalizability of underlying 
unspecified in other respects model. This results in a set of relations 
between low energy couplings of Z' to fermions that reduces in an 
essential way the number of parameters to be fitted in experiments. 
On this ground the observables which uniquely pick out the Z' boson 
in leptonic processes are introduced and the data of LEP experiments 
analyzed. The Z' couplings to leptons and quarks are estimated at 
95% confidence level. These estimates may serve as a guide for ex- 
periments at the Tevatron and/or LHC. A comparison with other 
approaches and results is given. 
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1 Introduction 



The precision test of the standard model (SM) at the LEP gave a possibihty 
not only to determine all the parameters and particle masses at the level of 
radiative corrections but also afforded an opportunity for searching for signals 
of new heavy particles beyond the energy scale of it. On the base of the LEP2 
experiments the low bounds on parameters of various models extending the 
SM have been estimated and the scale of new physics was obtained [H |2[ [3] . 
Although no new particles were discovered, a general believe is that the 
energy scale of new physics to be of order 1 TeV, that may serve as a guide 
for experiments at the Tevatron and LHC In this situation, any information 
about new heavy particles obtained on the base of the present day data is 
desirable and important. 

Numerous extended models include the Z' gauge boson - massive neutral 
vector particle associated with the extra f/(l) subgroup of an underlying 
group. Searching for this particle as a virtual state is widely discussed in 
the literature (see for references [H [S])- In the content of searching for Z' 
at the LHC and the ILC an essential information and prospects for future 
investigations are given in lectures [6]. Such aspects as the mass of Z', 
couplings to the SM particles, Z — Z' mixing and its influence in various 
processes and particles parameters, distinctions between different models are 
discussed in details. We shall turn to these papers in what follows. As 
concerned a searching for Z' in the LEP experiments and the experiments at 
Tevatron it was carried out mainly in a mo del- dependent way. A wide 
class of popular models has been investigated and low bounds on the mass 
mzi were estimated (see, [H El E])- As it is occurred, the low masses are 
varying in a wide energy interval 400-1800 GeV dependently on a specific 
model. These bounds are a little bit different in the LEP and Tevatron 
experiments. In this situation a model-independent analysis is of interest. 

In the papers [H [9l [TOl [11] of the present authors a new approach for the 
model-independent search for Z'-boson was proposed which, in contrast to 
other model-independent searches, gives a possibility to pick out uniquely 
this virtual state and determine its characteristics. The corresponding ob- 
servables have also been introduced and applied to analyze the LEP2 experi- 
ment data. Our consideration is based on two constituents: 1) The relations 
between the effective low-energy couplings derived from the renormalization 
group (RG) equation for fermion scattering amplitudes. We called them 
the RG relations. Due to these relations, a number of unknown Z' param- 
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eters entering the amplitudes of different scattering processes considerably 
decreases. 2) When these relations are accounted for, some kinematics prop- 
erties of the amplitude become uniquely correlated with this virtual state 
and the Z' signals exhibit themselves. 

The RG relations allow to introduce observables correlated uniquely with 
the Z'-boson. Comparing the mean values of the observables with the neces- 
sary specific values, one could arrive at a conclusion about the Z' existence. 
The confidence level (CL) of these values has been estimated and adduced 
in addition. Without taking into consideration the RG relations the deter- 
mination of Z'-boson requires a supplementary specification due to a larger 
number of different couplings contributing to the observables. Similar situ- 
ation takes place in the "hehcity model fits" of LEP Collaborations PQ El E] 
when different virtual states contribute to each of the specific models (AA, 
VV, and so on). Therefore these fits had the goal to discover any signals of 
new physics independently of the particular states which may cause devia- 
tions from the SM. Note that the LEP Collaborations saw no indications of 
new contact four fermion interactions in these fits. 

In Refs. [HI EDI [H] the one-parametric observables were introduced and 
the signals (hints in fact) of the Z' have been determined at the \o CL in 
the e^e~ /i"*"/!" process, and at the 2a CL in the Bhabha process. The 
Z' mass was estimated to be 1-1.2 TeV. An increase in statistics could make 
these signals more pronounced and there is a good chance to discover this 
particle at the LHC. 

In Ref. [12j the updated results of the one-parameter fit and the com- 
plete many-parametric fit of the LEP2 data were performed with the goal to 
estimate a possible signal of the Z'-boson with accounting for the final data 
of the LEP collaborations DELPHI and OPAL [2, 3J. Usually, in a many- 
parametric fit the uncertainty of the result increases drastically because of 
extra parameters. On the contrary, in our approach due to the RG relations 
between the low-energy couplings there are only 2-3 independent parameters 
for the investigated leptonic scattering processes. As it was showed in Ref. 
[12], an inevitable increase of confidence areas in the many-parametric space 
was compensated due to accounting for all accessible experimental informa- 
tion. Therefore, the uncertainty of the many-parametric fit was estimated as 
the comparable with previous one-parametric fits in Refs. [TOl [H]. In this 
approach the combined data fit for all lepton processes is also possible. 

From the results obtained on the searching for Abelian Z' within the 
LEP experiment data set we conclude that it is insufficient for convincing 
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discovery of this particle as the virtual state. In this situation it is reasonable 
to analyze the data by using the neural network approach which is able to 
make a realistic prognoses for the parameters of interest. This investigation 
was done within the two parametric global fit of the LEP2 data on the 
Bhabha scattering process. As the result of all these considerations we derive 
at the 2(7 CL the characteristics of the Z' (the vector v and axial-vector a 
couplings of the Z' with SM leptons and the Z — Z' mixing). The Z' mass 
is also estimated. Due to the universality of the a we also derived the model 
independent estimate of the Z' axial- vector couplings to quarks, Og = a. Note 
that the hints for the Z' have been determined in all the processes considered 
that increases the reliability of the signal. These results may serve as a good 
input into the future LHC and ILC experiments and used in various aspects. 
To underline the importance of them we mention that there are many tools 
at the LHC for the identification of Z' . But many of them are only applicable 
if Z' is relatively hght. The knowledge of the Z' couplings to SM fermions 
also have important consequences. 

The paper is organized as follows. In sect. 2 we give a necessary informa- 
tion about the description of Z' at low energies. In sects. 3-5 we discuss the 
origin of the RG relations, their explicit forms for the case of the heavy Z' and 
consequences of the relations for scattering processes investigated. In sect. 6 
the cross sections and the observables to pick out uniquely the virtual Z' in 
the e+e" fx^ ,t~^t~ processes are given. The fits of data are described 

and discussed. Then in sect. 7 the same is present for the Bhabha process 
e+e" —>■ e+e~. The one parametric and two parametric fits are discussed. In 
sect. 8 the analysis of this process is carried out by using the neuron network 
approach. The criteria for training the network are introduced which guaran- 
tee the 2(7 CL deviations of the data from the model containing the SM with 
extra Z'. The obtained parameters of the Z' practically coincide with that 
of derived in the one parameter analysis. In this way we determine the char- 
acteristics of the Z' coming from the LEP experiments. In sect. 9 we discuss 
the role of the present model-independent analysis for the LHC experiments. 
The discussion and comparison with results of other approaches are given in 
sect. 10. In the Appendix we describe the two-mass- scale Yukawa model and 
analyze in detail how the decoupling of the loop contributions due to heavy 
virtual states is realized when the mixing of fields is taken into consideration. 
This point is an essential element of the approach developed. 
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2 The Abelian Z' boson at low energies 



Let us adduce a necessary information about the Abelian Z'-boson. This 
particle is predicted by a number of grand unification models. Among them 
the and 5*0(10) based models [13] (for instance, LR, and so on) are 

often discussed in the literature. In all the models, the Abelian Z'-boson is 
described by a low-energy f/(l) gauge subgroup originated in some symmetry 
breaking pattern. 

At low energies, the Z'-boson can manifest itself by means of the couplings 
to the SM fermions and scalars as a virtual intermediate state. Moreover, 
the Z-boson couplings are also modified due to a Z-Z' mixing. In principle, 
arbitrary effective Z' interactions to the SM fields could be considered at 
low energies. However, the couplings of non-renormalizable types have to be 
suppressed by heavy mass scales because of decoupling. Therefore, significant 
signals beyond the SM can be inspired by the couplings of renormalizable 
types. Such couplings can be derived by adding new f/(l)-terms to the 
electroweak covariant derivatives D"^^ in the Lagrangian [HI |T2] (review, 

HE]) 

f ^.W^; - %B,Y,, - %B,Y,^ h (1) 




+ iY.^nl'\^y^-^9B,Qi-^B,Ys^^fn. 
Ir 

where summation over all the SM left-handed fermion doublets, leptons and 
quarks, fi = {fu)L, {fd)L, and the right-handed singlets, fn = {fu)R, (Jd)R, 
is understood. Qj denotes the charge of / in positron charge units, F/^ = 
diag(Yy„, F/^), and 1/^ = —1 for leptons and 1/3 for quarks. 

For general purposes we derive the RG relations for the Z' beyond the SM 
with two light Higgs doublets (THDM) [8]. Z' interactions with the scalar 
doublets can be parametrized in a model-independent way as follows. 



E 



(2) 



In these formulas, g, g', g are the charges associated with the SU{2) l, U (l)y, 
and the Z' gauge groups, respectively, are the Pauli matrices, Y^^ = 
diag(Y'0^ J, Y^^ 2) is the generator corresponding to the gauge group of the Z' 
boson, and Y^^ is the f/(l)y hypercharge. 
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The Yukawa Lagrangian can be written in the form 



fL i=l 

+GuAfLmfu)R + {fu)R<PrfL]) , (3) 

where (/)^ = zcr2(/)* is the charge conjugated scalar doublet. 

The Lagrangian ([2]) leads to the Z-Z' mixing. The Z-Z' mixing angle Oq 
is determined by the coupling as follows 

9o=~^^^^^f^^Y, + o(4-], (4) 

where 6w is the SM Weinberg angle, and aem is the electromagnetic fine 
structure constant. Although the mixing angle is a small quantity of or- 
der m^?, it contributes to the Z-boson exchange amplitude and cannot be 
neglected at the LEP energies. 

In what follows we will also use the Z' couplings to the vector and axial- 
vector fermion currents defined as 

- 9 ^ , «/ - 9 ^ • [0) 

The Lagrangian ([1]) leads to the following interactions between the fermions 
and the Z and Z' mass eigenstates: 

Czf-f = ^Zjr[iv^z+l'af^)oos9o+ 

,5 



+ {vf + -f a/) sin 6*0 /, 

r 

2' 



^z'f-f = \z'Jj>^Uvf + j'af)cos9o- 



^?^ + 7^4z)sin^o /, (6) 



where / is an arbitrary SM fermion state; f j^, are the SM couplings of 
the Z-boson. 

Since the Z' couplings enter the cross-section together with the inverse 
Z' mass, it is convenient to introduce the dimensionless couplings 



which can be constrained by experiments. 

Low energy parameters i, ^^^2? '^lj, ^rj must be fitted in exper- 
iments. In most investigations they were considered as independent ones. 
In a particular model, the couplings V^^ i, YJ/,. 2, Ylj, Yrj take some specific 
values. In case when the model is unknown, these parameters remain poten- 
tially arbitrary numbers. However, this is not the case if one assumes that 
the underlying extended model is a renormalizable one. In the papers [Bl E] 
it was shown that these parameters are correlated due to renormalizability. 
We called them the RG relations. Since this notion is a key-point of our 
consideration, we discuss it in detail. 



3 Renormalization group relations 

What is RG relation? Generally speaking, this is a correlation between low 
energy parameters of interactions of a heavy new particle with known light 
particles of the SM following from the requirement that full unknown yet 
theory extending SM is to be renormalizable. 

Strictly speaking, RG relations are the consequence of two constituencies: 

1. RG equation for a scattering amplitude; 

2. Decoupling theorem. 

The latter one describes the modification of both the RG operator 
d d d d 

and an amplitude at the energy threshold A of new physics. Here, and 
7x-functions correspond to all the charges \a and fields and masses X of the 
underlying theory. 

The RG equation for a scattering amplitude / reads, 

where / accounts for as intermediate states either the light or heavy virtual 
particles of the full theory. The standard usage of the RG equation is to 
improve the amplitude by solving this equation for the operator V calculated 



8 



in a given order of perturbation theory. However, to search for heavy virtual 
particles, we will use Eq. in another way. 

First we note that for any renormalizable theory, the RG equation is just 
identity, if / and T> are calculated in a given order of loop expansion. In this 
case Eq. ([9]) expresses the well known fact that the structure of the divergent 
term coincides with the structure of the corresponding term in a tree-level 
Lagrangian. 

For example, in massless QED, the tree-level plus one-loop one-particle- 
irreducible vertex function describing scattering of electron in an external 
electromagnetic field A, T = T^^^ + T^^\ is 





If we calculate the RG operator in one-loop order 

where I3^^\ 'yA\llp'^ S'^^ the beta-function and the anomalous dimensions of 
electromagnetic and electron fields, respectively, and apply it to F, we obtain 

^ r«=fe^4-27i^^-7i^^l r(°)+0(e^). (11) 



Then, accounting for the values of 

/?a) = ^ ^« = ^ ^« = ^ (12) 

and the factor e in F^°\ we observe that the first and the last terms in the 
r.h.s. cancel. Since /x-dependent term in F^^^ is F|^^) = ^f;^ log/x^, we see that 
Eq. fllll) is identity in the order O(e^). 

Next important point is that in a theory with different mass scales the 
decoupling of heavy-loop contributions at the threshold of heavy masses. A, 
results in the following property: the running of all functions is regulated by 
the loops of light particles. Therefore, the (3 and 7 functions at low energies 
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are determined by the SM particles, only. This fact is the consequence of the 
decouphng theorem (THl H?] . 

The decouphng results in the redefinition of parameters at the scale A 
and removing heavy-particle loop contributions from RG equation [181 E] ■ 

A A^ 

Aa = A, + aj^ Jog — + 6^ Jog2 — + ■ • • , (13) 

X ^ x(l + a,JogA + ,,Jog^^ + ...), 

where Aq and X denote the parameters of the SM. They are calculated as- 
suming that no heavy particles are excited inside loops. 

The matching between both sets of parameters Aa, X and Aa, X is chosen 
at the normalization point /i ~ A, 

Aa|^j=A = Aa|;x=A, -^aU=A = -^a|/x=A- (14) 

The differential operator V in the RG equation is in fact unique; the 
apparently different T) in both theories are the same! 

Note that if a theory with different mass scales is specified one can freely 
replace the parameters Aa, X and Aa, X by each other [8]. 

An example of the derivation and the main features of the RG relations 
are shown in the Appendix for a simple model with different mass scales. 

If underlying theory is not specified, the set of Aa, X is unknown. The low 
energy theory consists of the SM plus the effective Lagrangian generated by 
the interactions of light particles with virtual heavy particle states. The low 
energy parameters A^ of these interactions are arbitrary numbers which must 
be constrained by experiments. By calculating the RG operator T) and the 
scattering amplitudes of light particles in this 'external field' in a chosen order 
of loop expansion, it is possible to obtain the model-independent correlations 
between A^ . These are just the RE relations. 



4 The RG relations for Z' boson couplings 

Let us derive the correlations between Y^.^^, ^^^,2? ^lj, Yrj appearing due to 
renormalizability of the underlying theory containing Z' . 

In our case, the RG invariance of the vertex leads to the equation 
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where 

^ = wr- = 7^ + E/5.|- -E7x^Tr%, (16) 

dlog/i alog/i „ aAa ^ alogX 

and 

rfA rflogX 

= -r. , 7x = --7^ (17) 

a log /i a log /i 

are computed with taking into account the loops of light particles. 

Now we derive the RG relations following from the one-loop consideration. 
In accordance with the previous sections, the one-loop RG equation for the 
vertex function reads 



/^/^ + ^^« (frflf^) = 0, (18) 

o log /i mz' \ mz' J 

where F^^^' and Fj^' the tree-level and one-loop contributions to the 

fermion-Z' vertex. D*^^) = I]/?^^-*^ — l^x d\o x one-loop level part 

of the RG operator. 

To calculate these functions, only the divergent parts of the one-loop 
vertices are to be calculated. The corresponding diagrams are shown in 
Fig. [H The fermion anomalous dimensions 7^'* can be calculated by using 
the diagrams in Fig. [2l Then, Eq. lfTSl) leads to algebraic equations for the 
parameters F^^ i, y^, 2, Yl,/, and Yrj which have two sets of solutions [8]: 

Y^.,. = F^,, = -r^,, = -n, (19) 

Ylj + Ylj' = 0, Yrj = 0, 

and 

y<t>r.^ = Y^2.^=Y^a^Yt>^ (20) 

Ylj = Ylj*, Yaj = Y,j + 2T,fY^. 

Here / and /* are the partners of the SU{2)l fermion doublet (/* = ui, v* = 
I, g* = qd and = qu), T^f is the third component of weak isospin. 

The first of these relations describes the Z' boson analogous to the third 
component of the SU{2)l gauge field. The couplings to the right-handed 
singlet are absent. 

The second relation corresponds to the Abelian Z'. In this case the SM 
Lagrangian appears to be invariant with respect to the U{1) group associated 
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Figure 1: The diagrams contributing to the divergent parts of the Z'ff vertex 
at the one-loop leveL 
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Figure 2: The diagrams contributing to the fermion anomalous dimension at 
the one-loop level. 



with the Z' . The last relation in Eq. (l20l) ensures the Lvuk. Eq-© is to be 
invariant with respect to the f/(l) transformations. 

Introducing the Z' couplings to the vector and axial-vector fermion cur- 
rents ([5]), the last line in Eq. fl20l) yields 

Vf-af= Vf. - ttf* , af = T3fgY^. (21) 

The couplings of the Abelian Z' to the axial-vector fermion current have a 
universal absolute value proportional to the Z' coupling to the scalar doublet. 

These relations are model independent. In particular, they hold in all the 
known models containing the Abelian Z'. The most discussed models are 
derived from the Eg group (the so called LR, X'i' models). The tree-level Z' 
couplings to the SM fermions in the models are shown in Table 1. 

The Eg-symmetry breaking scheme 

Eg ^ SO(IO) X U(l)^, SO(IO) ^ SU(3)e x SU(2)z. x SU(2)^ x U(1)b-l. 

leads to the so called left-right (LR) model. Another scheme, 

Eg ^ SO(IO) X U(l)^, ^ SU(5) X U(l)x X U(l)^, 

predicts the Abelian Z', which is a linear combination of the neutral vector 
bosons ip and 

Z' = X cos P + ipsinP 
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Table 1: The Z' couplings to the SM fermions in the most discussed Eg-based 
models. 



/ 




LR 












V 


q COS /3 sin /3 


q COS (3 1 sin (3 


1 

2a 


1 

2« 


e 


cos /3 sin /3 


r) COS f3 

^ VTo 


a 
2 


1 « 
a 2 


Qu 


cos l3 sin f3 

VTo Ve 





a 
2 




Qd 


cos /3 sin 13 


r) COS /3 


a 
2 


1 « 
3a 2 



with the mixing angle /?. If we suppose only one Z' boson at low energies, 
the ip boson should be much heavier than the x field. In this case, the field 
ip is decoupled and /? — ^ 0. As it is seen, both the LR and the X'i^ models 
(with (3 = to avoid two Z' bosons with the same scale of masses) satisfy 
the RG relations (!20|) except for neutrinos. Let us explain this discrepancy. 
It is usually supposed in theories based on the Eg group that the Yukawa 
terms responsible for generation of the Dirac masses of neutrinos must be 
set to zero [TH]. Therefore, the terms proportional to the Yukawa couplings 
vanish in the renormalization group equation, and there are no RG relations 
for the Z' interactions with the neutrino axial-vector currents. In this case 
the couplings a^, given in Table 1 are not restricted by the RG relations. 

5 Implication of the RG relations 

LEP collaborations applied model dependent search for Z' and have obtained 
low bounds on the mass rriz' > 400—800 GeV dependently on a specific model 

In our analysis, we consider the SM with the additional effective Z' inter- 
actions ([1]), ([2]), ([3]) as a low energy theory. The parameters a/,f/ and mz' 
must be fitted in experiments. The RG relations give a possibility: 
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1. reduce the number of fitted parameters; 

2. determine kinematics of the processes; 

3. introduce observables which uniquely pick out the Z' signals. 

The RG relations (!20|) influences the Z — Z' mixing (jl]) . The axial- vector 
coupling determines also the coupling to the scalar doublet and, consequently, 
the mixing angle. As a result, the number of independent couplings is sig- 
nificantly reduced. 

In what follows, both types of the RG relations f|T9l) and fl20l) will be used 
in order to search for signals of the Z' gauge boson. 

6 Z' search in e^e~ ii^ ii~ ^t^t^ processes 
6.1 The differential cross section 

Let us consider the processes e+e^ l^l^ (/ = /i, r) with the non-polarized 
initial- and final-state fermions. In order to introduce the observable which 
selects the signal of the Abelian Z' boson we need to compute the differential 
cross-sections of the processes up to the one-loop level. 

The lower-order diagrams describe the neutral vector boson exchange in 
the s-channel {e~^e~ ^ V* ^ V = A,Z,Z'). As for the one- loop 

corrections, two classes of diagrams are taken into account. The first one 
includes the pure SM graphs (the mass operators, the vertex corrections, 
and the boxes) . The second set of the one-loop diagrams improves the Born- 
level Z'-exchange amplitude by "dressing" the Z' propagator and and the Z'~ 
fermion vertices. We assume that Z' states are not excited inside loops. Such 
an approximation means that the Z'-boson is completely decoupled. Then, 
the differential cross-section consists of the squared tree-level amplitude and 
the term from the interference of the tree-level and the one-loop amplitudes. 
To obtain an infrared-finite result, we also take into account the processes 
with the soft-photon emission in the initial and final states. 

In the lower order in m^? the Z' contributions to the differential cross- 
section of the process e+e^ l^l^ are expressed in terms of four-fermion 
contact couplings, only. If one takes into consideration the higher-order cor- 
rections in m^?, it becomes possible to estimate separately the Z'-induced 
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contact couplings and the Z' mass [T^. In the present analysis we keep the 
terms of order O(m^f) to fit both of these parameters. 

Expanding the differential cross-section in the inverse Z' mass and ne- 
glecting the terms of order 0[m~^\ we have 

% — 1 j — 1 

1 i 3 k 

+ E E E E Ci,fc„(s, z)x,x^x^,Xr,, (22) 

i=\ j=l k=\ n=l 

where the dimensionless quantities 

C = (xi,X2,a;3,X4,a;5,a;6,X7) = (a,We,^^/.,^^r,^^d,w^,^^6) (23) 

are introduced. Since the axial-vector couplings of the Abelian Z' boson are 
universal, we use the shorthand notation a = 0,^- In what follows the index 
I = fi,T denotes the final-state lepton. 

The coefficients A, B, C are determined by the SM couplings and masses. 
Each factor may include the tree-level contribution, the one-loop correction 
and the term describing the soft-photon emission. The factors A describe 
the leading-order contribution, whereas others correspond to the higher order 
corrections in m^?. 



6.2 The observable 

To take into consideration the correlations (jlj) we introduce the observable 
(Ti{z) defined as the difference of cross sections integrated in some ranges of 
the scattering angle 9 [9l [10] : 

aAz) = C -^dcosO - r -^dcosO, (24) 
^ ^ Jz dcosO J-i dcosO ' ^ ^ 

where z stands for the cosine of the boundary angle. The idea of introduc- 
ing the z-dependent observable is to choose the value of the kinematic 
parameter z in such a way that to pick up the characteristic features of the 
Abelian Z' signals. 
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The deviation of the observable from its SM value can be derived by the 
angular integration of the differential cross-section and has the form: 



7 i 



Aaiiz) 



i=l 0=1 
7 i j _ 

i=l j=l k=l n=l 



(25) 



Then let us introduce the quantity Act (z) = a {z) — asM (-z), which owing 
to the relations fl20l) can be written in the form 



aNf g% 



rrir 



F^{z,s)Y;^ + 2Fl{z,s)T,fYLjYL 



+ 2F/(z, s)T^fYLj% + Fi{z, s)YL,eY^ 



(26) 



The factor functions F/{z, s) depend on the fermion type through the \Qf\, 
only. In Fig. [3] they are shown as the functions of z for ^/s = 500 GeV. The 
leading contributions to F/{z, s), 



Fi{z,s) 



Fl{z,s) 



Fi(z,s) 



rrir. 



l-z^~\Qf\ {3z + 



+ 



ml 



- 1-^^ +- 3. + .^ (4|g;|-l) + 0^ 



rrir. 



Fi(z.s) 



\Qf\ [l-3z-z'-z']+0 



rrir. 



(27) 



are given by the Z' exchange diagram e~e+ — — //, since the Z-Z' 
mixing contribution to the Z exchange diagram is suppressed by the factor 
m|/s. 

From Eqs. (I27|) one can see that the leading contributions to the leptonic 
factors FI, F2, -F3 are found to be proportional to the same polynomial in 
z. This is the characteristic feature of the leptonic functions F- originating 
due to the kinematic properties of fermionic currents and the specific values 
of the SM leptonic charges. Therefore, it is possible to choose the value of 
z = z* which switches off three leptonic factors F^, F^, F^ simultaneously. 
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Figure 3: The leptonic functions Fq (the sohd curve), Fl (the long-dashed 
curve), F2 (the dashed curve), and Fg (the dotted curve) at ^/s = 500 GeV. 



Moreover, the quark function in the lower order is proportional to the 
leptonic one and therefore is switched off, too. As is seen from Fig. [3l the 
appropriate value of z* is about ~ 0.3. By choosing this value of z* one can 
simplify Eq. (126|) . It is also follows from Eq. (l26l) that neglecting the factors 
Fl, F2, Fl one obtains the sign definite quantity Aai{z*) ~ ~ a^. 

There is an interval of values of the boundary angle, at which the factors 
A[^, B{i, and C'nn at the sign-definite parameters a^, a?(, and contribute 
more than 95% of the observable value. It gives a possibility to construct the 
sign-definite observable Aai{z*) < by specifying the proper value of z*. 

In general, one could choose the boundary angle z* in different schemes. 
If just a few number of tree-level four-fermion contact couplings are consid- 
ered, one can specify z* in order to cancel the factor at the vector-vector 
coupling. However, if one-loop corrections are taken into account there are 
a large number of additional contact couplings. So, we have to define some 
quantitative criterion F{z) to estimate the contributions from sign-definite 
factors at a given value of the boundary angle z. Maximizing the criterion, 
one could derive the value z*, which corresponds to the sign-definite observ- 
able Aai{z*). Since the observable is linear in the coefficients A, B, and C, 
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Table 2: The boundary angle z* and the coefficients in the observable A(T/(z*) 
for the scattering into and r pairs at the one-loop level. 



v/i, GeV 






V" 






r 






z* 


An 


fill 


Ciiii 


z* 


All 


fill 


Ciiii 


130 


0.450 


-729 


-1792 


-19636 


0.460 


-687 


-1664 


-25782 


136 


0.439 


-709 


-1859 


-16880 


0.442 


-688 


-1779 


-20784 


161 


0.400 


-643 


-2183 


-6890 


0.400 


-625 


-2097 


-10993 


172 


0.390 


-619 


-4099 


-4099 


0.391 


-601 


-2263 


-8382 


183 


0.383 


-599 


-2545 


-1334 


0.385 


-571 


-2402 


-7580 


189 


0.380 


-586 


-2635 


-495 


0.380 


-568 


-2533 


-5135 


192 


0.380 


-579 


-2681 


-63 


0.380 


-562 


-2578 


-4769 


196 


0.380 


-571 


-2745 


-528 


0.379 


-554 


-2640 


-4272 


200 


0.378 


-564 


-2811 


-1137 


0.378 


-547 


-2704 


-3761 


202 


0.376 


-560 


-2845 


-1448 


0.377 


-543 


-2736 


-3501 


205 


0.374 


-555 


-2897 


-1923 


0.374 


-548 


-2834 


-1292 


207 


0.372 


-552 


-2932 


-2245 


0.372 


-544 


-2868 


-1010 



we introduce the following criterion, 



|y4n| + ^b\Bii\ + cuclCnnl 



all A 







Bij 




Cijkn 




all B 




all C 





(28) 



where the positive 'weights' ub ^ C and tuc ~ e take into account the order 
of each term in the inverse Z' mass. 

The numeric values of the 'weights' ujb and uc can be taken from the 
present day bounds on the contact couplings [Ij. As the computation shows, 
the value of z* with the accuracy 10~^ depends on the order of the 'weight' 
magnitudes, only. So, in what follows we take cub ~ 4 x 10"^ and uc ~ 
4 X 10-^. 

The function z*{s) is the decreasing function of the center-of-mass energy. 
It is tabulated for the LEP2 energies in Table O The corresponding values 
of the maximized function F are within the interval 0.95 < F < 0.96. 

Since A[^{s, z*) < 0, fin(s, z*) < and C[i^^{s, z*) < 0, the observable 



Aai{z*)= A[,{s 



Cfi 



11 



s, z 



111 



^s, z )a 



(29) 
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is negative with the accuracy 4-5%. Since this property follows from the 
RG relations (jlj) for the Abelian Z' boson, the observable A(T/(2*) selects the 
model-independent signal of this particle in the processes e^e~ — > It 
allows to use the data on scattering into \x\x and rr pairs in order to estimate 
the Abelian Z' coupling to the axial-vector lepton currents. 

Although the observable can be computed from the differential cross- 
sections directly, it is also possible to recalculate it from the total cross- 
sections and the forward-backward asymmetries. The recalculation proce- 
dure has the proper theoretical accuracy. Nevertheless, it allows to reduce 
the experimental errors on the observable, since the published data on the 
total cross-sections and the forward-backward asymmetries are more precise 
than the data on the differential cross-sections. 

The recalculation is based on the fact that the differential cross-section 
can be approximated with a good accuracy by the two-parametric polynomial 
in the cosine of the scattering angle z: 

^ = ^?!i^ + (i + ^^W + ^,, + *W. (30) 

az az 

where Si{z) measures the difference between the exact and the approximated 
cross-sections. The approximated cross-section reproduces the exact one in 
the limit of the massless initial- and final-state leptons and if one neglects 
the contributions of the box diagrams. 

Performing the angular integration, it is easy to obtain the expression for 
the observable: 

Aai{z*) = ai{z*) - af^(^*) = (1 - z*')vi - ^z*{3 + z*') + ~6i{^*), (31) 

y 

and for the total and the forward-backward cross-sections: 

A T T T,SM 8A , f / T \ 

Act, = ai -ai' = — + 6i{-l), 
Aai = Gi -ai =ni + di{0). (32) 
Then, the factors jSi and rji can be eliminated from the observable: 

Aa,{z*) = {l-z*')Aar-^z*i3 + z*')AaJ + ^i. (33) 
The quantity ^i, 

= Si{z*)-{l-z*')m + ^z*{3 + z*')~6i{-l), (34) 
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measures the theoretical accuracy of the approximation. 

The forward-backward cross-section is related to the total one and the 
forward-backward asymmetry by means of the following expression 

AaP = Aa,^ A^ + aP^ AA™. (35) 

As the computation shows, 6i{z*) ~ OMAai{z*), 6i{0) ^ 0.007Ao-f^, and 
6i{-l) ~ -O.OTAa^T at the LEP2 energies. Taking into account the exper- 
imental values of the total cross-sections and the forward-backward asym- 
metries at the LEP2 energies {Aaf ~ O.lpb, aj'^^ ~ 2.7pb, AAP =i 0.04, 
A™ ~ 0.5), one can estimate the theoretical error as C,i — O.OOSpb. At the 
same time, the corresponding statistical uncertainties on the observable are 
larger than O.OGpb. Thus, the proposed approximation is quite good and 
can be successfully used to obtain more accurate experimental values of the 
observable. 



6.3 Data fit 

To search for the model-independent signals of the Abelian Z'-boson we will 
analyze the introduced observable Aai{z*) on the base of the LEP2 data 
set. In the lower order in m^? the observable fl29|) depends on one flavor- 
independent parameter a^, 

Aa^\z*) = A[,{s, z*)a' + C[,,,{s, z*)a\ (36) 

which can be fitted from the experimental values of A(T^(z*) and Ao"^(z*). As 
we noted above, the sign of the fitted parameter (a^ > 0) is a characteristic 
feature of the Abelian Z' signal. 

In what follows we will apply the usual fit method based on the likelihood 
function. The central value of is obtained by the minimization of the x^- 
f unction: 

[Adf r^*) - AcT*^(2*)f 

= ? s^zM' ■ 

where the sum runs over the experimental points entering a data set chosen. 
The la CL interval (&i,62) for the fitted parameter is derived by means of 
the likelihood function £(a^) oc exp[— x^(a^)/2]. It is determined by the 
equations: 

/ ' C{e')de' = 0.68, C{bi) = £(62). (38) 

Jbi 
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To compare our results with those of Refs. PP we introduce the contact 
interaction scale 

= 4m|a-l (39) 

This normalization of contact couplings is admitted in Refs. . We use again 
the likelihood method to determine a one-sided lower limit on the scale A at 
the 95% CL. It is derived by the integration of the likelihood function over 
the physically allowed region a? > 0. The strict definition is 



A = 2mz(e* 



-1/2 



C{e')de' = 0.95 / C{e')de'. 



(40) 



We also introduce the probability of the Abelian Z' signal as the integral 
of the likelihood function over the positive values of d^: 



P 



L{e')de'. 



(41) 



Actually, the fitted value of the contact coupling originates mainly 
from the leading-order term in the inverse Z' mass in Eq. fl29p . The analysis 
of the higher-order terms allows to estimate the constraints on the Z' mass 
alone. Substituting in the observable fl29l) by its fitted central value, one 
obtains the expression 



^fitted 



4 

fitted ' 



(42) 



which depends on the parameter ( = m|/m|,. Then, the central value of 
this parameter and the corresponding la CL interval are derived in the same 
way as those for a?. 

To fit the parameters and ( we start with the LEP2 data on the total 
cross-sections and the forward-backward asymmetries [1]. The corresponding 
values of the observable Aai{z*) with their uncertainties 6ai{z*) are calcu- 
lated from the data by means of the following relations: 



Aai{z*) 



A™ (i 



^*2 



^(S + Z*' 
4 V 



+ l-z 



*2 



O", 



T,SM 
I 



AA 



FB 
I 5 



*\2 



Sai{z*) 



Ar (l - z*^) -'-{3 + 



T\2 



(San 



1 - z 



*2 



T,SM 
I 



FB\2 



(43) 
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Table 3: The contact coupling with the 68% CL uncertainty, the 95% CL 
lower limit on the scale A, the probability of the Z' signal, P, and the value 
of C = "ml/ml, as a result of the fit of the observable recalculated from the 
total cross-sections and the forward-backward asymmetries. 



Data set 




A, TeV 


P 


c 


jjLjJ, 
TT 

fifi and TT 


U.UUUUOOO_o 0000486 

U.UUUUZOO_o.ooo0639 
U.UUUUiOO_o.oo00387 


16.4 
17.4 
19.7 


0.77 
0.34 
0.63 


0.009 ±0.278 
-0.001 ±0.501 
0.017 ±0.609 



We perform the fits assuming several data sets, including the nn, tt, and 
the complete /i/x and rr data, respectively. The results are presented in Table 
[31 As is seen, the more precise /i/x data demonstrate the signal of about la 
level. It corresponds to the Abelian Z'-boson with the mass of order 1.2-1.5 
TeV if one assumes the value of a = ^^/47r to be in the interval 0.01-0.02. 
No signal is found by the analysis of the rr cross-sections. The combined fit 
of the /i/i and rr data leads to the signal below the la CL. 

Being governed by the next-to- leading contributions in m^?, the fitted 
values of ( are characterized by significant errors. The /i/i data set gives the 
central value which corresponds to m^/ ~ 1.1 TeV. 

We also perform a separate fit of the parameters based on the direct 
calculation of the observable from the differential cross-sections. The exper- 
imental uncertainties of the data on the differential cross-sections are of one 
order larger than the corresponding errors of the total cross-sections and the 
forward-backward asymmetries. These data also provide the larger values of 
the contact coupling a^. As for the more precise /i/i data, three of the LEP2 
Collaborations demonstrate positive values of a^. The combined is also 
positive and remains practically unchanged by the incorporation of the rr 
data. 

As it was mentioned in the previous section, the indirect computation of 
the observable from the total cross-sections and the forward-backward asym- 
metries inspires some insufficient theoretical uncertainty about 2% of the 
statistical one. It also increases the statistical error because of the recalcu- 
lation procedure. Nevertheless, the uncertainty of the fitted parameter 
within the recalculation scheme is of one order less than that for the direct 
computation from the differential cross-sections. This difference is explained 
by the different accuracy of the available experimental data on the differential 
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Figure 4: The factors at the Z' couphngs in the differential cross-section of 
the Bhabha process. 



and the total cross-sections. 



7 Search for Z' in e^e e+e process 
7.1 The differential cross-section 

In our analysis, as the SM values of the cross-sections we use the quantities 
calculated by the LEP2 collaborations [21 [3l [201 EI] • They account for either 
the one-loop radiative corrections or initial and final state radiation effects 
(together with the event selection rules, which are specific for each experi- 
ment). As it is reported by the DELPHI Collaboration, there is a theoretical 
error of the SM values of about 2%. In our analysis this error is added to the 
statistical and systematic ones for all the Collaborations. As it was checked, 
the fit results are practically insensitive to accounting for this error. 

The deviation from the SM is computed in the improved Born approx- 
imation. This approximation is sufficient for our analysis leading to the 
systematic error of the fit results less than 5-10 per cents. 

The deviation from the SM of the differential cross-section for the process 
e^t' V^t~ can be expressed through various quadratic combinations of 
couplings a = Oe, Ve, Vf^, Vr- For the Bhabha process it reads 



da da^^ , „„„, , , av. 



, , -/r(^)^ + /r(^)^ + /rw^, (44) 

dz dz rriyi my, rriy. 



where the factors are known functions of the center-of-mass energy and the 
cosine of the electron scattering angle z plotted in Fig. HI The deviation of 
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the cross-section for e~^e 



fi~^^ (r+r ) processes has a similar form 



da da"^^ 




dz dz 



+ 



(45) 



Eqs. (I44l) - (l45l) are our definition of the Z' signal. 

Note again that the cross-sections in Eqs. ( l44l) - (H5l) account for the re- 
lations (l62il through the functions fi{z), fs^z), f4{z), since the coupling 
(the mixing angle 9q) is substituted by the axial coupling constant a. Usu- 
ally, when a four-fermion effective Lagrangian is applied to describe physics 
beyond the SM [22j, this dependence on the scalar field coupling is neglected 
at all. However, in our case, when we are interested in searching for signals 
of the Z'-boson on the base of the effective low-energy Lagrangian ([I])-®, 
these contributions to the cross-section are essential. 

7.2 One-parameter fit 

The factor /|^(-z) is positive monotonic function of z (see Fig. 4 for the center- 
of-mass energies ^/s = 200 GeV. The same behavior is observed for higher 
energies). Such a property allows one to choose f2^{z) as a normalization 
factor for the differential cross section. Then the normalized deviation of the 
differential cross-section reads [TT] 



and the normalized factors are shown in Fig 5. Now these factors are finite at 
z ^ 1. Each of them in a special way influences the differential cross-section. 

1. The factor at is just the unity. Hence, the four-fermion contact cou- 
pling between vector currents, w^, determines the level of the deviation 
from the SM value. 

2. The factor at a? depends on the scattering angle in a non-trivial way. It 
allows to recognize the Abelian Z' boson, if the experimental accuracy 
is sufficient. 



da 

dz 



m| da 



infl^z) dz 



(46) 
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Figure 5: Factors Fa{\/s, z) (solid) and Fav{\fs-i z) (dashed) in the normahzed 
deviation of the differential cross-section da/dz for ^/s — 200 and 500 GeV. 
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3. The factor at av results in small corrections. 

Thus, effectively, the obtained normalized differential cross-section is a 
two-parametric function. In the next sections we introduce the observables 
to fit separately each of these parameters. 

7.3 Observables to pick out v'^ 

To recognize the signal of the Abelian Z' boson by analyzing the Bhabha pro- 
cess the differential cross-section deviation from the SM predictions should 
be measured with a good accuracy. At present, no such deviations have 
been detected at more than the la CL. In this situation it is resonable to 
introduce integrated observables allowing to pick out Z' signals by using the 
most effective treating of available data. The observables should be sensitive 
to the separate Z' couplings. This admits of searching for the Z' signals in 
different processes as well as to perform global fits. 

The normalized deviation of the differential cross-section f l46|) is (effec- 
tively) the function of two parameters, and v"^. We are going to introduce 
the integrated observables which determine separately the four-fermion cou- 
plings and v"^ [Tl]. 

Let us first proceed with the observable for v"^. After normalization the 
factor at the vector-vector four-fermion coupling becomes the unity. Whereas 
the factor at is a sign- varying function of the cosine of the scattering angle. 
As it follows from Fig. 5, for the center-of-mass energy 200 GeV it is small 
over the backward scattering angles. So, to measure the value of the 
normalized deviation of the differential cross-section has to be integrated 
over the backward angles. For the center-of-mass energy 500 GeV the factor 
at is already a non- vanishing quantity for the backward scattering angles. 
The curves corresponding to intermediate energies are distributed in between 
two these curves. Since they are sign- varying ones at each energy point some 
interval of z can be chosen to make the integral to be zero. Thus, to measure 
the Z' coupling to the electron vector current v'^ we introduce the integrated 
cross-section (146|) 



where at each energy the most effective interval [zq, Zq + A2;] is determined 
by the following requirements: 




(47) 
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1. The relative contribution of the couphng v"^ is maximal. Equivalently, 
the contribution of the factor at a? is suppressed. 

2. The length /S.z of the interval is maximal. This condition ensures that 
the largest number of bins is taken into consideration. 

The relative contribution of the factor at v"^ is defined as 



Az + 



-Az 



Fadz 



+ 



(48) 



and shown in Fig. 6 as the function of the left boundary of the angle interval, 
zq, and the interval length, Az. In each plot the dark area corresponds to 
the observables which values are determined by the vector- vector coupling v'^ 
with the accuracy > 95%. The area reflects the correlation of the width of 
the integration interval Az with the choice of the initial zq following from the 
mentioned requirements. Within this area we choose the observable which 
includes the largest number of bins (largest A^;). The corresponding values 
of Zq and Az are marked by the white dot on the plots in Fig. 6. As the 
carried out analysis showed, the point Zq is shifted to the right with increase 
in energy whereas Az remains approximately the same. 

From the plots it follows that the most efficient intervals are 

-0.6 <z < 0.2, = 200 GeV, 

-0.3 <z< 0.7, ^/s = 500 GeV. (49) 

Therefore the observable (H7I) allows to measure the Z' coupling to the elec- 
tron vector current v"^ with the efficiency > 95%. 

Fitting the LEP2 final data with the one-parameter observable, we find 
the values of the Z' coupling to the electron vector current together with 
their la uncertainties: 



ALEPH 


v! = 


-0.11 ±6.53 X 10~ 


DELPHI 


vl = 


1.60 ± 1.46 X 10" 


-4 


L3 


vl = 


5.42 ±3.72 X 10" 


-4 


OPAL 


vl = 


2.42 ± 1.27 X 10" 


-4 


Combined 


t = 


2.24 ±0.92 X 10" 


-4 



As one can see, the most precise data of DELPHI and OPAL collaborations 
are resulted in the Abelian Z' hints at one and two standard deviation level, 
correspondingly. The combined value shows the 2o" hint, which corresponds 
to 0.006 < \vA < 0.020. 
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Figure 6: Relative contribution of the factor at v"^ to the observable ay as 
the function of the left boundary of the angle interval, zq^ and the interval 
length, A^;, at the center-of-mass energy 200 and 500 GeV. The shaded areas 
correspond to the contributions > 95% (dark) , from 90% to 95% (midtone) , 
and < 90% (hght). 
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7.4 Observables to pick out 

In order to pick the axial- vector coupling one needs to eliminate the domi- 
nant contribution coming from . Since the factor at in the dd/dz equals 
unity, this can be done by summing up equal number of bins with positive and 
negative weights. In particular, the forward-backward normalized deviation 
of the differential cross-section appears to be sensitive mainly to a^, 

/•^max da /"O , da 
(Tfb = / dz- / dz — 

Jo dz J-zm^^ dz 

= Fa^FBO^ + Fav, FBdV. (50) 

The value ^max is determined by the number of bins included and, in fact, 
depends on the data set considered. The LEP2 experiment accepted e^e~ 
events with \z\ < 0.72. In what follows we take the angular cut z^ai^x = 0.7 
for definiteness. 

The efficiency of the observable is determined as: 





Fa,FB 




\Fa,FB 


1 + 





(51) 



It can be estimated as k = 0.9028 for the center-of-mass energy 200 GeV and 
K = 0.9587 for 500 GeV. Thus, the observable 

apB = 0.224a2 - 0.024at;, ^/^ = 200 GeV, 

aFB = 0.472a2 - 0.020av, ^ = 500 GeV (52) 

is mainly sensitive to the Z' coupling to the axial- vector current a^. 

Consider a usual situation when experiment is not able to recognize the 
angular dependence of the differential cross-section deviation from its SM 
value with the proper accuracy because of loss of statistics. Nevertheless, a 
unique signal of the Abelian Z' boson can be determined. For this purpose the 
observables J^°'^^^{da/dz)dz and a-pB must be measured. Actually, they are 
derived from the normalized deviation of the differential cross-section. If the 
deviation is inspired by the Abelian Z' boson both the observables are to be 
positive quantities simultaneously. This feature serves as the distinguishable 
signal of the Abehan Z' virtual state in the Bhabha process for the LEP2 
energies as well as for the energies of future electron-positron collider ILC 
(> 500 GeV). The observables fix the unknown low energy vector and axial- 
vector Z' couplings to the electron current. Their values have to be correlated 
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with the bounds on and v"^ derived by means of independent fits for other 
scattering processes. 

We estimated the observable fl52l) related to the value of a^. Since in the 
Bhabha process the effects of the axial-vector coupling are suppressed with 
respect to those of the vector coupling, we expect much larger experimental 
uncertainties for a^. Indeed, the LEP2 data lead to the huge errors for of 
order 10~^ — 10~^ The mean values are negative numbers which are too large 
to be interpreted as a manifestation of some heavy virtual state beyond the 
energy scale of the SM. 

Thus, the LEP2 data constrain the value of v"^ at the 2cr CL which could 
correspond to the Abelian Z' boson with the mass of the order 1 TeV. In 
contrast, the value of is a large negative number with a significant exper- 
imental uncertainty. This can not be interpreted as a manifestation of some 
heavy virtual state beyond the energy scale of the SM. 



7.5 Many-parameter fits 



As the basic observable to fit the LEP2 experiment data on the Bhabha 
process we propose the differential cross-section 



^^Bhabha ^^Bhabha,S'Af 



dz 



dz 



(53) 



where i runs over the bins at various center-of-mass energies ^/s. The fi- 
nal differential cross-sections measured by the ALEPH (130-183 GeV, |20j). 
DELPHI (189-207 GeV, [5]), L3 (183-189 GeV, [21j), and OPAL (130-207 
GeV, [2j) collaborations are taken into consideration (299 bins). 

As the observables for e~^e~ — >■ n'^fi~,T'^T~ processes, we consider the 
total cross-section and the forward-backward asymmetry 



e+i-,SM 



A 



FB 



A 



e+i-,SM 

FB 



(54) 



where i runs over 12 center-of-mass energies y/s from 130 to 207 GeV. We 
consider the combined LEP2 data [Ij for these observables (24 data entries 
for each process). These data are more precise as the corresponding differ- 
ential cross-sections. Our analysis is based on the fact that the kinemat- 
ics of s-channel processes is rather simple and the differential cross-section 
is effectively a two-parametric function of the scattering angle. The total 
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cross-section and the forward-backward asymmetry incorporate complete in- 
formation about the kinematics of the process and therefore are an adequate 
alternative for the differential cross-sections. 

The data are analysed by means of the fit [12j. Denoting the observ- 
ables (I53l) - (l5^ by 0"^, one can construct the x^-function, 

2 

(55) 

where cx^^ and 5a are the experimental values and the uncertainties of the 
observables, and a^^ are their theoretical expressions presented in Eqs. f H4l) - 
(l45|l . The sum in Eq. (155!) refers to either the data for one specific process 
or the combined data for several processes. By minimizing the x^-function, 
the maximal-likelihood estimate for the Z' couplings can be derived. The x^- 
function is also used to plot the confidence area in the space of parameters 
a, Vei and Vr- Note that in this way of experimental data treating all 
the possible correlations are neglected. We believe that at the present stage 
of investigation this is reasonable, because the Collaborations have never 
reported on this possibility. 

For all the considered processes, the theoretic predictions af^ are linear 
combinations of products of two Z' couplings 

j:Q,A„ (56) 

{a^, vl, aVe, VeVf,, VeVr, aV^, ttV^}, 

where Cij are known numbers. In what follows we use the matrix notation 
cr*^ = af^, a^^ = erf"", C = Cij, A = Aj. The uncertainties 6ai can be 
substituted by a covariance matrix D. The diagonal elements of D are ex- 
perimental errors squared. Da = (Sa^^)'^, whereas the non-diagonal elements 
are responsible for the possible correlations of observables. The x^-function 
can be rewritten as 

X^{A) = (a^"-(T*^)TD-i(a^"-(T*^) 

= (a""-CA)^D-^((T""-CA), (57) 

where the upperscript T denotes the matrix transposition. 
The x^-function has a minimum, Xmin; 

A = {C^D~^C)''C^D~^a'''' (58) 
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err 



5ai 



a. 



th 



A, 



corresponding to the maximum-likelihood values of Z' couplings. From Eqs. 
dSZ]), dSHD we obtain 



x'(^)-xLn = {A-Afb-\A-A\ 

B = {C^D-^C)-\ (59) 

Usually, the experimental values a'^^ are normal-distributed quantities 
with the mean values o"*^ and the covariance matrix D. The quantities A, 
being the superposition of a^"^, also have the same distribution. It is easy to 
show that A has the mean values A and the covariance matrix D. 

The inverse matrix is symmetric and can be diagonalized. The num- 
ber of non-zero eigenvalues is determined by the rank (denoted M) of D~^. 
The rank M equals to the number of linear-independent terms in the observ- 
ables cr*^. So, the right-hand-side of Eq. fl59l) is a quantity distributed as 
with M degrees of freedom (d.o.f.). Since this random value is independent 
of A, the confidence area in the parameter space (a, Ve, v^, Vr) corresponding 
to the probability [3 can be defined as [23] : 

x'<xLn + XcL,;3W, (60) 

where Xcl,/3(^) is the /3-level of the x^-distribution with M d.o.f. 

In the Bhabha process, the Z' effects are determined by 3 linear-independent 
contributions coming from a?, v'^, and aVe (M = 3). As for the e~^e~ — >■ 
r+r~ processes, the observables depend on 4 linear- independent terms 
for each process: a^, VeVf^, VeO,, dv^ for e~^e~ — > and a?, VeVr, Vga, dVr 

for e^e~ t^t~ (M = 4). Note that some terms in the observables for 
different processes are the same. Therefore, the number of d.o.f. in the 
combined fits is less than the sum of d.o.f. for separate processes. Hence, 
the predictive power of the larger set of data is not drastically spoiled by 
the increased number of d.o.f. In fact, combining the data of the Bhabha 
and e"'"e~ fi~^fi~ (r+r^) processes together we have to treat 5 linear- 
independent terms. The complete data set for all the lepton processes is 
ruled by 7 d.o.f. As a consequence, the combination of the data for all the 
lepton processes is possible. 

The parametric space of couplings (a, Ve, v^, Vr) is four-dimensional. 
However, for the Bhabha process it is reduced to the plane (a, Ve), and to 
the three-dimensional volumes (a, Ve, v^), (a, Ve, Vr) for the e+e" 
and e"^e" r^r^ processes, correspondingly. The predictive power of data 
is distributed not uniformly over the parameters. The parameters d and 
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-0.02 0.02 a 

Figure 7: The 95% CL areas in the (a, Ve) plane for the Bhabha, e+e^ — >■ 
//"•"/i", and e+e~ t^t^ processes. 

Ve are present in all the considered processes and appear to be significantly 
constrained. The couplings or Vr enter when the processes e+e^ fJ'^f^~ 
or e~^e~ t^t~ are accounted for. So, in these processes, we also study the 
projection of the confidence onto the plane (a,fe). 

The origin of the parametric space, a = Ve = 0, corresponds to the absence 
of the Z' signal. This is the SM value of the observables. This point could 
occur inside or outside of the confidence area at a fixed CL. When it lays 
out of the confidence area, this means the distinct signal of the Abelian Z'. 
Then the signal probability can be defined as the probability that the data 
agree with the Abelian Z' boson existence and exclude the SM value. This 
probability corresponds to the most stringent CL (the largest Xcl) which 
the point a = Ve = is excluded. If the SM value is inside the confidence area, 
the Z' boson is indistinguishable from the SM. In this case, upper bounds on 
the Z' couplings can be determined. 

The 95% CL areas in the (a, Ve) plane for the separate processes are plot- 
ted in Fig. [71 As it is seen, the Bhabha process constrains both the axial- 
vector and vector couplings. As for the e'^e' — » yU^/x^ and e"'"e^ t^t^ pro- 
cesses, the axial-vector coupling is significantly constrained, only. The confi- 
dence areas include the SM point at the meaningful CLs, so the experiment 
could not pick out clearly the Abelian Z' signal from the SM. An important 
conclusion from these plots is that the experiment significantly constrains 
only the couplings entering sign-definite terms in the cross-sections. 

The combination of all the lepton processes is presented in Fig. [HI There is 
no visible signal beyond the SM. The couplings to the vector and axial- vector 
electron currents are constrained by the many-parameter fit as l^el < 0.013, 
\a\ < 0.019 at the 95% CL. If the charge corresponding to the Z' interactions 
is assumed to be of order of the electromagnetic one, then the Z' mass should 
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Figure 8: The projection of the 95% CL area onto the {a,Ve) plane for the 
combination of the Bhabha, e~^e~ — fJ^^n" , and e^e^ t^t^ processes. 

be greater than 0.67 TeV. For the charge of order of the SM SU{2) l couphng 
constant m^/ > 1.4 TeV. One can see that the constraint is not too severe to 
exclude the Z' searches at the LHC 

Let us compare the obtained results with the one-parameter fits. As one 
can see, the most precise data of DELPHI and OPAL collaborations are 
resulted in the Abelian Z' hints at one and two standard deviation level, 
correspondingly. The combined value shows the 2a hint, which corresponds 
to 0.006 < l^el < 0.020. On the other hand, our many-parameter fit con- 
strains the Z' coupling to the electron vector current as l^el < 0.013 with no 
evident signal. Why does the one-parameter fit of the Bhabha process show 
the 2(7 CL hint whereas there is no signal in the two-parameter one? Our 
one-parameter observable accounts mainly for the backward bins. This is in 
accordance with the kinematic features of the process: the backward bins 
depend mainly on the vector coupling v1, whereas the contributions of other 
couplings are kinematically suppressed (see Fig. 4). Therefore, the difference 
of the results can be inspired by the data sets used. To clarify this point, we 
perform the many-parameter fit with the 113 backward bins {z < 0), only. 
The minimum, x^j^ = 93.0, is found in the non-zero point \a\ = 0.0005, 
Ve = 0.015. This value of the Z' coupling Ve is in an excellent agreement 
with the mean value obtained in the one-parameter fit. The 68% confidence 
area in the (a, Vg) plane is plotted in Fig. O There is a visible hint of the 
Abelian Z' boson. The zero point d = Vg = (the absence of the Z' boson) 
corresponds to = 97.7. It is covered by the confidence area with 1.3a 
CL. Thus, the backward bins show the 1.3a hint of the Abelian Z' boson in 
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Figure 9: The 68% CL area in the (a, Ve) plane from the backward bins of the 
Bhabha process in the LEP2 experiments (the shaded area). The hatched 
area is the 68% CL area from the LEP 1 data on the Bhabha process. 

the many-parameter fit. So, the many-parameter fit is less precise than the 
analysis of the one-parameter observables. 

At LEPl experiments [24j the Z-boson couplings to the vector and axial- 
vector lepton currents {gy, Qa) were precisely measured. The Bhabha process 
shows the la deviation from the SM values for Higgs boson masses mn > 114 
GeV (see Fig. 7.3 of Ref. [24j). This deviation could be considered as the 
effect of the Z-Z' mixing. It is interesting to estimate the bounds on the Z' 
couplings following from these experiments. 

Due to the RG relations, the Z-Z' mixing angle is completely determined 
by the axial-vector coupling a. So, the deviations of gy, qa from their SM 
values are governed by the couplings a and Ve, 

9v - g!^ = -m.06ave, gA - gT = 49.06a^ (61) 

Let us assume that the total deviation of theory from experiments follows due 
to the Z-Z' mixing. This gives an upper bound on the Z' couplings. In this 
way one can estimate whether the Z' boson is excluded by the experiments 
or not. 

The la CL area for the Bhabha process from Ref. [21] is converted into 
the (a, Ve) plane in Fig. [91 The SM values of the couplings correspond to the 
top quark mass mj = 178 GeV and the Higgs scalar mass ttih = 114 GeV. 
As it is seen, the LEPl data on the Bhabha process is compatible with the 
Abelian Z' existence at the la CL. The axial- vector coupling is constrained 
as \d\ < 0.005. This bound corresponds to < 2.5 x 10~^, which agrees with 
the one-parameter fits of the LEP2 data for e'^e~ — ^ fj,'^fi~,T^T~ processes 
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Table 4: The summary of the fits of the LEP data for the dimensionless 
contact couphngs ([7]). 



Data 




a2 


LEPl 


e-e+, 68% CL 




(1.25 ± 1.25) X 10-5 


LEP2, one-parameter fits 


e-e+, 68% CL 


(2.24 ±0.92) X 10-4 




/i/i, 68% CL 




(3.66;t86) X 10-5 


liH,TT, 68% CL 




(1.33ti|?) X 10-5 


LEP2, many-parameter fits 


e e+, /i/i, rr, 95% CL 


< 1.69 X 10-4 


< 3.61 X 10-4 


e~e~^ backward, 68% CL 


(2.25i^:o^?) X 10-4 


< 9.49 X 10-4 



{a? = 1.3 ± 3.89 X 10-5 at 68% CL). On the other hand, the vector couphng 
constant Ve is practically unconstrained by the LEPl experiments. 

For the convenience, in Table 4 we collect the summary of the fits of the 
LEP data in terms of dimensionless contact couplings (JTj). From the analysis 
carried out we come to conclusion that, in principle, the LEP experiments 
were able to detect the Z'-boson signals if the statistics had been sufficient. 

8 Z' hints within neural network analysis 

Since the actual LEP2 data set is not too large to detect Z' boson, one needs 
in the estimate of its parameters which could be used in future experiments. 
To determine them in a maximally full way we address to the analysis based 
on the predictions of the neural networks (for applications in high energy 
physics see, for example, P5|). The main idea of this approach is to constrain 
a given data set in such a way that an amount of it is considered as an 
inessential background and omitted. The remaining data are expected to 
give a more precise fit of the parameters of interest. 

We take into consideration the complete set of the differential cross sec- 
tions for the Bhabha process accumulated by all the LEP Collaborations and 
apply the following criteria to restrict the data [26] : 

1. As the signal we use the differential cross sections for the Bhabha pro- 
cess accounting for the SM plus Z' and calculated at 0.25 x 10-4 < 



37 



vl<Ax 10"^ and 0.25 x 10"^ <al<Ax lO"! 

Such a choice of parameters is motivated by the results obtained in the 
previous sections. The cross sections due to the Z' exchange diagrams 
were calculated with the RG relations been taken into consideration. 

2. As the background we use the deviations of the experimental differential 
cross sections from calculated for the SM plus Z' ones which are larger 
than redoubled uncertainties of LEP2 experimental data. 

The network trained with these criteria omits the events which correspond 
to the large deviations from the theoretical cross sections but accounts for 
the peculiarities proper the Z' existence. To construct and train the neural 
network we used the program MLPFit ^7]. The results of the carried out 
analysis based on the two parametric fit discussed in the previous sections 
demonstrate the 2cr CL hint for the Z' . For the vector coupling the neural net- 
work at the 2a CL predicts = 2.4±1 .99 x 10"^ [26] that is in agreement with 
the discussed above result derived in the one parametric fit. The obtained 
values of v1 correspond to the value of the mass mz' = 0.53 — 1.05 TeV, if the 
coupling g is of the order of the SM gauge couplings, g'^/ (47r) ~ 0.01 — 0.03. 
Thus, the carried out analysis demonstrates the hint of Z' boson which can 
be not too heavy. We conclude once again that the data set of the LEP 
experiments is not sufficient to detect the pronounced signal of this virtual 
particle. 



9 Search for Chiral in Bhabha process 

Let us turn to the analysis of the Bhabha process with the aim to search for 
the Chiral Z' gauge boson The Chiral Z' interacts with the SM doublets 
only that can be described by one parameter for each doublet (Yfi and F^). 
It is characterized by the constraints (fT9|) 



YLj = -YfL(Ts, Yrj = 0, Y^^ = -Y^a, (62) 

where is the Pauli matrix. 

Remind that in the Bhabha process it is convenient to use the normalized 
cross-sections fH6l): 



da rn?y da 7^-2 - - 7^-7 7^- 

'dz ^ '^kJ^^z) dZ ^ ""'^ ^ FavaeVe + A,</,^^e0 + Fa^^a^ 
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Since the Chiral Z' boson docs not interact with the right-handed species, the 
normahzed deviation of the differential cross-section from its SM prediction 
is determined by two factors, and Fj,,^, 

da mi ^ da ^ ^ t ^ 
where we define the dimensionless constants 

h = -JT — ^ = ~7r — s^'t'- 

The normahzation gives us two benefits. First, the obtained factors F{^/s, z) 
are finite for all values of the scattering angle z. Second, the experimental 
uncertainties for different bins become equalized that provides the statistical 
equivalence of different bins. The latter is important for the construction of 
integrated cross- sect ions. 

The Z-Z' mixing angle is determined by as follows, 

. mw sin 9w 7 
do ^ (/>, 

where aem is the fine structure constant. 



9.1 One-parameter fit for Chiral Z' 

The Chiral Z' boson does not interact with the right-handed species. The 
normalized deviation of the differential cross-section from its SM prediction, 

da/dz = FJI + FL4,le4>, 
is determined by two finite factors, and F^^, which are shown in Fig. 10. 

As it is seen, the four-fermion contact coupling contributes mainly 
to the forward scattering angles, whereas the Z-Z' mixing term affects the 
backward angles. At the LEP energies they can be of the same order of 
magnitude. The contribution of the mixing vanishes with the energy growth. 
We have a possibility to derive the effective experimental constraints on them 
without any additional restrictions. 

First, let us construct a one-parametric observable which is most preferred 
by the statistical treatment of data. As is clear, it is impossible to separate 
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Figure 10: The normalized factors and Fi^^ describing the Chiral Z' 
effects in the Bhabha process at \fs = 200 GeV (sohd hnes) and y/s = 500 
GeV (dashed hnes). 

the couphngs Zg and le(f) in any observable which is an integrated cross-section 
over some interval of z. However, the mixing contribution can be eliminated 
in the cross-section of the form (which is inspired by the forward-backward 
asymmetry) 



where the boundary value z* should be chosen to suppress the coefficient 
at le4>- The maximal value of the scattering angle z^^^^ is determined by a 
particular experiment. In this way we introduce the one-parametric sign- 
definite observable sensitive to J^. 

The LEP Collaborations DELPHI and L3 measured the differential cross- 
sections with Zmax = 0.72 [3l [21]. The set of boundary angles z* as well as 
the theoretic and experimental values of the observable are collected in Table 
5. The other LEP Collaborations - ALEPH and OPAL - used z^^^ = 0.9 
[2], [20]. The corresponding data are presented in Table 6. 

The standard x^"fif gives the following constraints for the coupling at 



ALEPH 


F 


= -0.00304 ±0.00176 


DELPHI 


''e 


= -0.00054 ±0.00086 


L3 


''e 


= 0.00109 ±0.00184 


OPAL 


''e 


= 0.00051 ± 0.00064 


Combined 


''e 


= -0.00004 ± 0.00048 




the 68% CL: 
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Tabic 5: The boundary angles z* and the theoretic and experimental values 
of the observable Act (2;*) at Zmax = 0.72. 



Vs, GeV 


z* 




AcT^^(;2*), DELPHI 


A(7^'^(;2*), L3 


183 


-0.245 


1742 I'i 




-2.38 ±7.03 


189 


-0.252 


1775 ll 


-4.28 ±3.36 


3.05 ± 3.65 


192 


-0.255 


1788 ll 


4.57 ±7.32 




196 


-0.259 


1806 ll 


3.77 ±4.33 




200 


-0.263 


1823 ll 


-1.95 ±4.05 




202 


-0.265 


1831 Tl 


1.31 ±5.54 




205 


-0.267 


1843 ll 


-4.09 ±3.89 




207 


-0.269 


1851 ll 


0.40 ± 3.33 





Table 6: The boundary angles z* and the theoretic and experimental values 
of the observable Acr(;2*) at ^^max = 0.9. 



yi, GeV 


z* 


A(7(Z*) 


A(j^^(z*), ALEPH 


A(7^^(z*), OPAL 


130 


-0.217 


2017 ll 


-12.40 ± 19.24 


-4.13 ±29.29 


136 


-0.266 


2092 ll 


-50.21 ± 16.64 


-34.18 ±31.58 


161 


-0.370 


2311 ll 


-15.90 ± 13.24 


-14.02 ±22.32 


172 


-0.400 


2398 ll 


-12.11 ± 12.50 


13.71 ± 17.84 


183 


-0.424 


2474 ll 


-1.51 ± 5.18 


11.04 ± 5.57 


189 


-0.435 


2512 ll 




-0.63 ± 3.28 


192 


-0.441 


2531 [2 




-3.48 ±9.85 


196 


-0.447 


2554 ll 




2.96 ± 5.09 


200 


-0.454 


2577 ll 




0.35 ± 4.68 


202 


-0.457 


2587 ll 




-2.87 ± 9.00 


205 


-0.461 


2604 ll 




5.88 ± 4.67 


207 


-0.464 


2614 ll 




-1.42 ± 3.46 
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Figure 11: The 95% CL area in Zg ~ plane. The final data of ALEPH 
130-183 GeV, DELPHI GeV 189-207, L3 183-189 GeV, and OPAL 130-207 
GeV are combined. 

Hence it is seen that the most precise data of DELPHI and OPAL collabora- 
tions give no signal of the Chiral Z' at the la CL. The combined value also 
shows no signal at the la CL. From the combined fit the 95% CL bound on 
the value of T1 can be derived, fg < 9 x 10""^. Supposing the Z' coupling con- 
stant g to be of the order of the electroweak one, g ~ 0.6, the corresponding 
Z' mass has to be larger than 0.5 TeV. 

9.2 Two parametric fit for Chiral Z' 

A complete two parametric fit of experimental data based directly on the 
differential cross-sections has been carried out in Ref. [28]. Due to only two 
independent couplings this fit is efficient. In the fitting the available final 
data for the differential cross-sections of the Bhabha process were used. The 
data set consists of 299 bins including the data of ALEPH at 130-183 GeV, 
DELPHI at 189-207 GeV, L3 at 183-189 GeV, and OPAL at 130-207 GeV 
[21 [21 [201 [21] • The fitting procedure is similar to that of discussed above for 
the Abelian Z' . The results can be summarized as follows. 

The parameter space of the Chiral Z' is the plane (Zg, 0)- The minimum 
of the x^-function, x^^^ = 237.29, is reached at zero value of Tg (~ 10""^) and 
almost independent of the value of (j) (the maximal-likelihood values of the 
couplings). The 95% CL area (xcl — 5.99) is shown in Figure 3. 




42 



As one can see, the zero point, /e = = (the absence of the Chiral Z' 
boson) is inside the confidence area. The value of in this point (238.62) 
is indistinguishable from the Xmin- other words, the set of experimental 
data cannot determine the signal of the Chiral Z'-boson. 

As also is seen, the value of le is constrained as le < 0.02 at the 95% CL. 
This upper bound is in an agreement with the corresponding result of the 
one-parameter fit (/g < 0.03). Thus, the Z' mass has to be larger than 0.75 
TeV, if the Z' coupling constant g is again supposed to be of the order of the 
electroweak one, g ~ 0.6. 

The fit of the differential cross-sections leads to a better accuracy for le 
than the fit of the integrated cross-sections based on the same data. Without 
accounting for the model-independent relations between the Z' couplings it 
is impossible to obtain such results. 

10 Model independent results and search for 
Z' at the LHC 

In this section we discuss all the assumptions giving a possibility to pick 
out the Z' signal and determine its characteristics in a model independent 
way. We also note the role of the present results for the future LHC and 
ILC experiments. As it was already stressed, in searching for this particle 
at the LEP and Tevatron a model dependent analysis was mainly used. The 
motivation for this was the different number of chiral fermions involved in 
different models (see, for example, [5]). In this approach, the low bounds 
on rriz' have been estimated and a smallness of the Z — Z' mixing was also 
observed. 

On the contrary, in our model independent approach the RG relations 
between the parameters of the effective low energy Lagrangians have been 
accounted for that gave a possibility to determine not only the bounds but 
also the mass and other parameters of the Z'. 

First we note all the assumptions used in our considerations. We ana- 
lyzed the four-fermion scattering amplitudes of order ~ m~^? generated by 
the Z' virtual states. The vertices linear in Z' were included into the effective 
low-energy Lagrangian. We also impose a number of natural conditions. The 
interactions of a renormalizable type are dominant at low energies ~ mw- 
The non-renormalizable interactions generated at high energies due to ra- 
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diation corrections are suppressed by the inverse heavy mass ~ and 
neglected. We also assumed that the SU{2)l x f/(l)y gauge group of the SM 
is a subgroup of the GUT group. As a consequence, all structure constants 
connecting two SM gauge bosons with Z' are to be zero. Hence, the interac- 
tions of gauge fields of the types Z'W^W~ , Z'ZZ, and other are absent at 
the tree level. Our effective Lagrangian is also consistent with the absence of 
the tree-level flavor-changing neutral currents (FCNCs) in the fermion sec- 
tor. The renormalizable interactions of fermions and scalars are described 
by the Yukawa Lagrangian ([3]) that accounts for different possibilities of the 
Yukawa sector without the tree-level FCNCs. These assumptions are quite 
general and satisfled in a wide class of Eq inspired models. 

Within these constraints for the low energy effective Lagrangian the RG 
correlations have been derived. Correspondingly, the model independent 
estimates of the mass m,z' and other parameters are regulated by the noted 
requirements. Therefore, the extended underlying model has also to accept 
them. 

In this regard, let us discuss the role of the obtained estimates for the 
LHC. As it is well known (see, for example, OE]), there are many tools at 
the LHC for Z' identification. But many of them are only applicable if this 
particle is relatively light. Our results are in favor to this case. 

Next important point is the determination of Z' couplings to the various 
SM fermions. As we have shown, the axial-vector couplings of the Z' to the 
SM fermions are universal and proportional to its coupling to the Higgs field. 
Hence we have obtained an estimate of the a = a/ couplings for both leptons 
and quarks. This is an essential input because experimental analyses for 
the LHC have mainly concentrated on being able to distinguish models and 
not on actual couplings. The vector coupling Ve was also estimated that, in 
particular, may help to distinguish the decay of the Z' resonance state to e~^e~ 
pairs. Since the couplings Oe and Ve were estimated there is a possibility to 
distinguish this process from the decay of the KK system. In the literature [U 
|29] on searching for the Z' it is also mentioned that the determination of the 
Z' couplings to fermions could be fulfilled channel by channel, Og, f e, Ve^b, cie,b, 
.... In almost all these considerations the relations between the parameters 
have not been taken into account. But this is very essential for treating of 
experimental data and introducing relevant observables to measure. 

Other parameter is the Z — Z' mixing, which is responsible for differ- 
ent decay processes and the effective interaction vertices generated at the 
LHC [HE]- It is also determined by the axial- vector coupling (see Eq. (^I^) 
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and estimated in a model independent way. Remind that in our analysis 
(in contrast to the approaches of the LEP Collaborations) the mixing was 
systematically taken into account. Its value is of the same order of magni- 
tude as the parameters that were fitted in experiments. Note also that the 
existence of other heavy particles with masses mx > ttlz' does not influence 
the RG relations which are the consequences of the necessary condition for 
renormalizability. In fact, this condition (the structure of a divergence gen- 
erated by radiation corrections coincides with that of the tree-level vertex) 
holds for each renormalizable type interaction. 

An important role of the model independent results for searching for Z' 
at the LHC and ILC consists, in particular, in possibility to determine the 
particle as a virtual state due to a large amount of relevant events. We men- 
tioned already that, in principle, LEP experiments were able to determine it 
if the statistics was sufficiently large. Experiments at the ILC will increase 
numerously the data set of interest. In fact, the observables, introduced in 
sects. 6 and 7 for picking out uniquely aj and couplings in the leptonic 
scattering process, are also effective at energies ^/s > 500 GeV and could be 
applied in future experiments at ILC. 

Other model independent methods of searching for the Z' as a resonance 
state are proposed in the literature (see Refs. [291 SOI EI]). We do not 
discuss them here because they take no relations between the parameters 
into consideration. Besides, the main goal of the present paper is to adduce 
model independent information about the Z' followed from experiments at 
low energies. Different aspects of Z' physics at the LHC are out of the scope 
of it. 

11 Discussion 

In this section we collect in a convenient form all the results obtained and 
make a comparison with other investigations on searching for Z' at low en- 
ergies. In fact, this is a large area to discuss. References to numerous results 
obtained in either model dependent or model independent approaches can 
be found in the surveys E]- Further subdivision can be done into the 
considerations accounting for any type correlations between the parameters 
of the low energy effective interactions and that of assuming complete inde- 
pendence of them. Because of a large amount of fitting parameters the latter 
are less predictable. 
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Table 7: The summary of the fits of the LEP data for the maximum likelihood 
values of the Z' couplings fl63|) to the SM fermions and of the Z — Z' mixing 
angle 9q. M = denotes the unknown value of the Z' mass in TeV units. 



Data 


l^ol, 


Kl, xio-i 


\a'f\, xlQ-^ 


A^, xlO-^ 


LEPl 


e e+ 


3.17M-1 




1.38M 


0.437 




LEP2, 


one-parameter fits 




e~e+ 




5.83M 








5.42M-1 




2.36M 


1.278 




3.27M-1 




1.42M 


0.464 


LEP2, many-parameter fits 


e e+, 2; < 




5.84M 







Now, for a convenience of readers let us present the results of fits of the 
Z' parameters in terms of the popular notations [ll|5]. The Lagrangian reads 



zff 



^Z'ff 



,SM 



+ A;)-7^(af + a;)I /, 



(63) 



with the SM values of the Z couplings 



,SM 



e(T3f -2Qf sin^ 6w 



3/ 



sin 9w cos 9 



w 



sin 9w cos 9 



w 



where e is the positron charge, Qj is the fermion charge in the units of e, 
T3J = 1/2 for the neutrinos and u-type quarks, and T-^f = —1/2 for the 
charged leptons and d-tjpe quarks. 

The results of the fits of the Z' couplings to the SM leptons obtained from 
the analysis of the LEP experiments are adduced in the Tables and Fig. 
[T^ Remind that due to the universality of the axial- vector coupling af the 
same estimates take also place for quarks. First of all, one parameter fits of 
LEP experiments as well as the many-parameter fit for the e~^e~ backward 
bins show the hints of the Z' boson at the l-2cr CL. Due to this fact, the fits 
allow to determine the maximum likelihood values of Z' parameters. In spite 
of uncertainties, these values can be used as a guiding line for the estimation 
of possible Z' effects in the LHC experiments. The maximum likelihood 
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Table 8: The summary of the fits of the LEP data for the confidence intervals 
for the Z' couplings fl63l) to the SM fermions and for the Z — Z' mixing angle 
6q. M = JZ?' denotes the unknown value of the Z' mass in TeV units. 



Data 


CL 


l^ol, 


Kl, xio-i 


\a'f\, xlQ-^ 




LEPl 


e e+ 


68% 


(0;4.48)M-i 


(0;1.95)M 


(0; 0.873) 


LEP2, one-parameter fits 


e e+ 


95% 




(2.46;7.87)M 








95% 


(0;10.39)M-^ 




(0;4.52)M 


(0; 4.694) 




95% 


(0;8.64)M-i 




(0;3.75)M 


(0; 3.244) 


LEP2, many-parameter fits 


e e^, 


95% 


(0;17.03)M-i 


(0;5.06)M 


(0; 7.40)M 


(0; 12.607) 


z <0 


68% 


(0;27.61)M-i 


(1.68;7.83)M 


(0; 12.00)M 


(0; 33.1288) 



values are given in Table [7l As is seen, different fits and processes lead to 
the comparable values of the Z' parameters. 

In Table E] we present the confidence intervals for the fitted parameters. 
With this Table one is able to estimate the uncertainty of the Z' couplings 
as well as the lower bounds on the parameters. 

In Fig. [12] the maximum likelihood values and the CL intervals are drawn 
for the different values of the Z' mass. All the plots exploit the same color 
scheme. The values excluded at 95% CL by the many-parameter fit of all the 
LEP2 leptonic processes e^e^ —>■ l^l^ are shown in gray. The 95% confidence 
intervals from the one-parameter fit of LEP2 e^e~ — ^ fJ'~^fJ'~ , t^t~ are drawn 
in pink with the maximum likelihood values as the dashed red line. The 
corresponding results with taking into account the /i"*"//" process only are 
shown in yellow with solid red line for the maximum likelihood values. The 
maximum likelihood values from the LEPl experiments are represented as 
dotted blue line. The 95% confidence interval from the one-parameter fit 
of the LEP2 Bhabha scattering is shown as the blue crosshatched region 
with the maximum likelihood values as the dashed blue line. Finally, the 
68% confidence interval and the maximum likelihood values from the many- 
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Figure 12: The maximum likelihood values and the confidence intervals for 
the Z — Z' mixing angle (^0)5 the vector coupling to electron current (fg), and 
the axial-vector couplings to the SM fermions (a'j) by the LEP 2 data. The 
values excluded at 95% CL by the many-parameter fit of e+e~ are 
shown in gray. The results of fits based on the one-parameter observables are 
shown in pink for e^e~ — > [i^ [i~ ^t^t" (with the maximum likelihood value 
as the dashed red line), in yellow for e^e~ — >■ [i^ [i~ (with the maximum 
likelihood value as the solid red line), and in blue for e^e~ — *■ e^e~ (with the 
maximum likelihood value as the dashed blue line). The maximum likelihood 
values and the Icr CL area for the many-parameter fit of backward bins of 
e^e~ e^e~ are also shown in green. The blue dotted lines correspond to 
the maximum likelihood values obtained from the LEPl data. 
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parameter fit of tlie backward bins of LEP2 Bhabha scattering are shown in 
green. 

Now we compare the above results with the ones of other fits accounting 
for the Z' presence. As it was mentioned in Introduction, LEP collaborations 
have determined the model dependent low bounds on the Z' mass which 
vary in a wide energy interval dependently on a model. The same has also 
been done for Tevatron experiments. The modern low bound is m^/ > 850 
GeV. It is also well known that though almost all the present day data 
are described by the SM [II [21 [Sj [2l] , the overall fit to the standard model 
is not very good. In Ref. [7] it was showed that the large difference in 
sin^ ^g'jj'* from the forward-backward asymmetry Ajj, of the bottom quarks 
and the measurements from the SLAG SLD experiment can be explained for 
physically reasonable Higgs boson mass if one allows for one or more extra 
[/(I) fields, that is Z' . A specific model to describe Z' physics of interest 
was proposed which introduces two type couplings to the hyper charge Y 
and to the baryon-minus-lepton number B — L. Within this model by using 
a number of precision measurements from LEPl, LEP2, SLD and Tevatron 
experiments the parameters ay and ob-l of the model were fitted. The 
presence of Z' was not excluded at 68% CL. The value of ay was estimated 
to be of the same order of magnitude as in our analysis and is comparable with 
values of other parameters detected in the LEP experiments. The erroneous 
claim that ay is two order less then the value derived from our Table 4 is, 
probably, a consequence of some missed factors. The upper limit on the mass 
was also obtained mz' < 2.6 TeV . 

These two analyzes are different but complementary. A common feature 
of them is an accounting for the Z' gauge boson as a necessary element of 
the data fits. The results are in accordance at 68-95% CL with the existence 
of the not heavy Z' which has a good chance to be discovered at Tevatron 
and/or LHC 

Appendix. RG relations in a theory with dif- 
ferent mass scales 

In this Appendix we are going to investigate the Yukawa model with a heavy 
scalar field x ^"^^ ^ light scalar field [32]. The goals of our investigation 
are two fold: 1) to derive the one- loop RG relation for the four-fermion scat- 
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Figure 13: Tree level contribution to the four-fermion amplitude. 



tering amplitude in the decoupling region and 2) to find out the possibility 
of reducing this relation in the equation for vertex describing the scattering 
of light particles on the external field when the mixing between heavy and 
light virtual states takes place. 

The Lagrangian of the model reads 

+p^\' + ij {td,^, -M-G^^- G^x) ^, (64) 

where i(j is a. Dirac spinor field, and A ^ m, M. 

Consider the four-fermion scattering tjjtp — ^ tpip. The S-matrix element 
at the one-loop level is given by 

i f dpi dpi 




IPR - <^^i<^<l,2 ^ Tc. ^2 \!,. Zl2 



m 



(t>2 ' 



Xlp {pi) [1 + r (p3, -p4 - Ps)] (ps) 

Xi> {pi) [1 + r (p2, -Pl - P2)\ ^ iP2) , (65) 

where s = [pi + ^2)^, -S'lPR is the contribution from the one-particle reducible 
diagrams shown in Figs. [TBlfT^ and Sbox is the contribution from box dia- 
grams. The one- loop polarization operator of scalar fields n^^^^ and the 
one-loop vertex function F are related to the Green functions as 



/ \ ^4^14*2 TT / \ "i 

/>102 \^) ~ ~ ~~2 "I 2 0102 \^) ~ 2 ' 



m^,, s - m^,, s -m. 



G^4>^{p,q) = -J2G^,D^^^(^q^^ S^{p) {1 + T{p,q))S^{-p -q), {66) 
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Figure 14: The one- loop level contribution to the one-particle reducible four- 
fermion amplitude. 



where is the spinor propagator in the momentum representation. 



Renormalization of the model 

The renormalized fields, masses and charges are defined as follows 



X 



-1/2 / V^O 

I Xo 



G 



Ml - 5M^, 



(67) 



where subscript marks the corresponding bare quantities. 

Using the dimensional regularization (the dimension of the momentum 
space is D = 4 — £) and the MS renormalization scheme one can compute 
the renormalization constants 



~ ~ 16^ ('^^ ^ ' 



8rin? 



87r% 



1 - 



87r% 



ioG 



G\ 



1 

1 - 



(gI + 6A) - 6GlM^ - pA' 
{gI + 6^) A'~QGlM^-pm^ 



(Gi+Gi)] [zrf. 



(68) 
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From Eq. 
loop level: 

7a 

7^ 



we obtain the appropriate /3 and 7 functions at the one- 



duG 

dlogK 

dlogK 
d log 



m 



q2 

-m? ^ r 



^ 5G, + 3Gp + 4 ^2_^2 W ' 



1 



c/loj 


2; K, 


47r2 


(ilog 




1 


c/loj 




47r2 


(ilog 


^ . 


= 1 ( 


(iloe 




327r2 V 



^^2m2-6M2 

^ + 6A - p— 



G: 



.A' 



- 6M2 ^ m2 



(69) 



Then, the S'-matrix element can be expressed in terms of the renormalized 
quantities. The contribution from the one-particle reducible diagrams be- 
comes 



'S'lPR — ^ G^jG 



+ 



XtPipl) \l + r(*^°) {p2, -Pi - P2)l ^(P2), 



(70) 



where the functions ^Ip^J^ T'^^'^^ are the expressions Il4>i(j,2 and F without 
the terms proportional to 1/e. Since the quantity ^box is finite, the renor- 
malization leaves it without changes. 

Introducing the RG operator at the one-loop level [H] 



V 



d 



d 



d 



d log K d log K d log K 

d 



d 



d 



In 



7a- 



d 



^ 'dG^ "' d\ogm? 91ogA2 d\og%l) 
we determine that the following relation holds for the S'-matrix element 



(71) 



(S'lPR + Sbox) 



^'^IPR I -T)(l) C(0) _ n 

a log K 
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(72) 



where the Sip^ and the S\p^ are the contributions to the Sipr at the tree 
level and at the one-loop level, respectively: 

six = + ^1 X (73) 



d log K 47r2 



Gl 



•£ |_ X 

m? s — 



Gl (pA2 - 6Am2 + Gl {6M^ - s 

(s — m^)^ 
2GlGli6M'-s) 



{s — m?) (s — A^) 
Gl [pm^ - e^A' + Gl {QM^ - s] 
(s-A2)2 



(74) 



The first term in Eq. ([71]) is originated from the one-loop correction to 
the fermion-scalar vertex. The rest terms are connected with the polarization 
operator of scalars. The third term describes the one-loop mixing between 
the scalar fields. It is canceled in the RG relation (1721) by the mass-dependent 
terms in the f3 functions produced by the non-diagonal elements in Z^. 

Eq. (1721) is the consequence of the renormalizability of the model. It 
insures the leading logarithm terms of the one-loop S'-matrix element to re- 
produce the appropriate tree-level structure. In contrast to the familiar treat- 
ment we are not going to improve scattering amplitudes by solving Eq. (!72|) . 
We will use it as an algebraic identity implemented in the renormalizable 
theory. Naturally if one knows the explicit couplings expressed in terms of 
the basic set of parameters of the model, this RG relation is trivially fulfilled. 
But the situation changes when the couplings are represented by unknown 
arbitrary parameters. In this case the RG relations are the algebraic equa- 
tions dependent on these parameters and appropriate (3 and 7 functions. In 
the presence of a symmetry the number of (3 and 7 functions is less than the 
number of RG relations. So, one has non trivial system of equations relating 
the unknown couplings. For example, such a scenario is realized for the gauge 
coupling. Although the considered simple model has no gauge couplings, we 
are able to demonstrate the general procedure of deriving the RG relations. 
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Decoupling of the heavy field 

At energies s <^ the heavy scalar field x is decoupled. This means, 
that the four-fermion scattering amplitude is described by the model with 
no heavy field x plus terms of the order s/A^. At the tree level, this is the 
obvious consequence of the expansion of the heavy scalar propagator 



1 1 



A2 A2 



01^ 



(75) 



which is resulted in the effective contact four-fermion interaction in Eq. (1731) 

Ccs = -aipi/j X ijji/j, a = -^. (76) 
So, the tree level contribution to the scattering amplitude becomes 



c.(0) 
"-•iPR 



s — rn? \A^ 



ip^ X ^ip, (77) 



and the lowest order effects of the heavy scalar in the decoupling region are 
described by the parameter a, only. 

Decoupling of heavy particles is present also at the level of radiative cor- 
rections. The radiative corrections are generally described by various loop 
integrals in the momentum space (the Passarino-Veltman functions). Con- 
sidering a Passarino-Veltman function with at least one heavy mass Lambda 
inside loop in the low-energy limit, one can see the following asymptotic 
behavior: the function splits into 1) possible energy-independent divergent 
part (including also log A) and 2) energy- dependent finite part which can be 
expanded by inverse powers of A and vanishes at the small energies. The 
important property is that the logA-term in the divergent part reproduces 
the logarithm of the cut-off scale. So, such a potentially large term has to be 
automatically absorbed by the renormalization at low energies and leads to 
no observable effects. However, if the renormalization is actually performed 
at high energies (as in the MS renormalization scheme) the potentially large 
log A-terms should be re-summed manually by the redefinition of the physical 
couplings and masses at the scale A. 

What is the form of the RG relations in the limit of large A? The method 
of constructing the RG equation in the decoupling region was proposed in 
|18j . It introduces the redefinition of the parameters of the model allowing 
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to remove all the heavy particle loop contributions to Eq. fl74l) . Let us define 
a new set of fields, charges and masses ip, G^, G-^, A, m, M 

01 - Gjfl.glog^^... 



57r^ A^ 



where dots stand for the higher powers of log A responsible for the decoupling 
at higher loop orders. 

The differential operator ( ffll) ban be rewritten in terms of these new 
low-energy parameters: 

dlogK dlogK y ^dG^ 

9 ^ d ^ d . . 

where /3 and 7 functions are obtained from the one-loop relations (169!) and 



- 1^ ( ~ ^^2 ^^^l 



' 47r2 I >^ I A2 i ^A2, 
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7* = (80) 

Hence, one immediately notices that (3 and 7 functions contain only the 
light particle loop contributions, and all the heavy particle loop terms are 
completely removed from them. The 5'-matrix element expressed in terms of 
new parameters satisfies the following RG relation 

V (5iPK + 5box) = Ifi^^ + ^W^Sj, = 0, (81) 
log ^L 

^Sr = I - « + O (i) ) X ^V^, (82) 



d log K Att"^ \ s — m? 

Gl (-QXw? + Gl (6M2 _ s 




s — m? 



(83) 



where a = G JlSr is the redefined effective four-fermion coupling. As one 



can see, Eq. f l83l) includes all the terms of Eq. fl7^ except for the heavy 
particle loop contributions. It depends on the low energy quantities -j/^, G<^, 
a. A, m, M. The first and the second terms in Eq. (1551) are just the one-loop 
amplitude calculated within the model with no heavy particles. The third 
and the fourth terms describe the light particle loop correction to the effective 
four-fermion coupling and the mixing of heavy and light virtual fields. 



Elimination of the one-loop scalar field mixing 

Due to the mixing term it is impossible to split the RG relation f lHTj) for 
the 5'-matrix element into the one for vertices. Hence, we are not able to 
consider Eq. flHTl) in the framework of the scattering of light particles on an 
external field induced by the heavy virtual scalar. But this is an important 
step in deriving the RG relation for EL parameters. Fortunately, there is a 
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simple procedure allowing to avoid the mixing in Eq. 0831) . The way is to 
incorporate the diagonalization of the leading logarithm terms of the scalar 
polarization operator into the redefinition of the (p, G^, G^. 



X 



C 



1/2 / V 

V X 



G^f 
G. 



1 

1 H — log — 

327r2 ^ A2 



-1/2 



^ G, 
\ G. 
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1/2 



GipGx 



87r2 (A2 - 
The appropriate $ functions 
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A2 







A2 - 6M2 




5GI 
167r2 



167r2 



2(jr^ + 3G-^G^p 



(85) 



contain no terms connected with mixing between light and heavy scalars. So, 
the fourth term in Eq. (1551) is removed, and the RG relation for the S'-matrix 
element becomes 



box j 



^'^iPR I -n(i) c(o) _ n 
a log K 



5 



(0) 
IPR 



s — 



- a + O 



(86) 
(87) 



agjp^R ^ Uy<U ( Gl 
d log K Atc"^ \ s — fh? 



G% (-6Am2 + G2 (6M2_^^^ 




At (5 = Eq. fl86|) is just the RG identity for the scattering amplitude 
calculated in the absence of the heavy particles. The terms of order a describe 
the RG relation for the effective low-energy four-fermion interaction in the 
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decoupling region. The last one can be reduced in the RG relation for the 
vertex describing the scattering of the light particle (fermion) on the external 
field a/5 substituting the virtual heavy scalar: 

V (Vdii)^) = -^V^i>ij + p(^) (v^V^^/i) = 0, (89) 

where 

"^dGp ^"Sloga ^""Slogm^ ^"^^log^' 
^^ = -^{^Gl + Oia)). (90) 

Eqs. fl5B]) - flM?]) is the main result of our investigation. One can derive them 
with only the knowledge about the low-energy couplings of heavy particle 
(!76|) and the Lagrangian of the model with no heavy particles. One also has 
to ignore all the heavy particle loop contributions to the RG relation and 
the one-loop mixing between the heavy and the light fields. Eqs. ( !86l) - (!90ll 
depend on the effective low-energy parameters, only. But as the difference 
between the original set of parameters and the low-energy one is of one-loop 
order, one may freely substitute them in Eqs. fl86l) - fl89|) . It is also possible 
to reduce the RG relation for scattering amplitudes to the one for vertex 
describing the scattering of light particles on the 'external' field induced by 
the heavy virtual particle. In fact, this result is independent on the specific 
features of the considered model. 
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